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' Abstract. The luminosity distribution of Galactic radio pulsars is believed to be log-normal in 
0^ . form. Applying this functional form to populations of pulsars in globular clusters, we employ 
' Bayesian methods to explore constraints on the mean and standard deviation of the function, as 
well as the total number of pulsars in the cluster. Our analysis is based on an observed number of 
pulsars down to some limiting flux density, measurements of flux densities of individual pulsars, 
as well as diffuse emission from the direction of the cluster. We apply our analysis to Terzan 
5 and demonstrate, under reasonable assumptions, that the number of potentially observable 
• pulsars is in a 95.45% credible interval of 133lgg^. Beaming considerations would increase the 
C^' true population size by approximately a factor of two. 
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. 1. Introduction 

K>" I Globular clusters have high core stellar number densities that favour the formation of 

■ lovif-mass X-ray binaries (LMXBs) that are believed to be the progenitors of miUisecond 
' pulsars (MSPs; Alpar et al. 1982). MSPs can be considered long-lived tracers of LMXBs, 

ly-j ■ so constraints on the MSP content provide unique insights into binary evolution and the 
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integrated dynamical history of globular clusters, vifhile determining the radio luminosity 
function of these pulsars helps shed light on their emission mechanism. 
. Faucher-Giguere & Kaspi (2006) have shown that the luminosity distribution of non- 

I ■ recycled Galactic pulsars appears to be log-normal in form. More recently, Bagchi et al. 
^ [ (2011) have verified that the observed luminosities of recycled pulsars in globular clusters 
are consistent with this result. Assuming, therefore, that there is no significant difference 
between the nature of Galactic and cluster populations, we use Bayesian techniques to 
^ ■ investigate some of the consequences that occur when one applies this functional form to 
■ ■ ' populations of pulsars in individual clusters. We are interested in the situation where we 
observe n pulsars with luminosities above some limiting luminosity. There is a family of 
luminosity function parameters (fx, a) and population sizes (N) that is consistent with 
this observation, and here we analyze the posterior probabilities of different members of 
this family given the data. In our case, the data are the individual pulsar flux densities 
that we call {Si}, the observed number of pulsars, n and the total diffuse flux density of 
the cluster, S'obs- 
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2. Bayesian parameter estimation 

Luminosity and flux density are related by the standard pseudo-luminosity equation 
L = S , where r is the distance to the pulsar (see Lorimer & Kramer 2005). This implies 
that the luminosity function is corrupted by uncertainties in distance. To mitigate this, 
we decided to perform our analysis initially in terms of the measured flux densities, and 
use a model of distance uncertainty to convert our results to the luminosity domain. We 
take the distance to all pulsars in a cluster to be the same. The log-normal in luminosity 
can then alternatively be written in terms of flux density. The probability of detecting a 
pulsar with flux density S in the range log S to log S + dilog S) is given by a log-normal 
in S as 

-I _ (log S-MS)^ 

p(log S) d(log S) = j=e ^'i d{\og S), (2.1) 

where S is in mJy, and fis and as are the mean and standard deviation of the flux density 
distribution. The probability of observing a pulsar above the limit S'min is then 

Po^s^f^ pilogS) dilogS) ^ \ erfc ( ^"^^J^"^^ ■ (2-2) 

First, we consider as data the measured flux densities of pulsars in the cluster, {Si}. 
Ideally, the survey sensitivity limit S'min can be taken as another datum, but its exact 
value is not always known, so we decided to parametrize S'min- The likelihood of observing 
a set of pulsars with fluxes {Si} is represented as 

" " 1 (log S.--f.s)^ 

nP.(l0gS,|MS,as,Smin)=n (2-3) 

,^1 ,^1 PobsfT5V27r 

where n is the number of observed pulsars in the cluster, and pobs is as given in Equa- 
tion (j2.2l) . Uncertainties in the flux density measurements are not considered here, but 
it has to be noted that it will have the effect of underestimating the credible intervals on 
our posteriors. 

To infer the total number of pulsars in the cluster, we follow Boyles et al. (2011) to 
take as likelihood the probability of observing n pulsars in a cluster with N pulsars, given 
by the binomial distribution 

p{n\N, AiS, OS, Smin) = ^,^^1^^^, P^bs (1 " Pobs)^""- (2.4) 

Next, we incorporate information about the observed diffuse flux from the direction of 
the cluster. We assume that all radio emission is due to the pulsars in the cluster, both 
resolved and unresolved. For the likelihood of measuring the diffuse flux S'obs, we choose 
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p(Sob.|iV,/iS,as) = j=e , (2.5) 

CTdiffV 27r 

where S'dift is the expectation of the total diffuse flux of a cluster whose flux density 
distribution is a log- normal with parameters /is and erg, and having N pulsars, and 
tTdifi is the standard deviation. Here, Sdig = N{S) and cTdiff = ^fN SD(S) where the 
expectation of S is given by (S) ioMs+i<T|in(io) ^Yie standard deviation of S, 
SD(S') = IQi's+h'^l \'!^m^iQ'^%Mm _ I -y^g Jo j^o^ consider the uncertainty in the 
diffuse flux measurement. The total likelihood, p(log S'i, n, Sobsl-^, /^s, erg, S'min) is the 
product of the three likelihoods computed above. 

The flux density distribution of pulsars in a cluster is not suitable for comparing 
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the populations in different clusters, as it depends on the distance to the cluster. So 
we transform the total likelihood obtained in the previous subsection to the luminosity 
domain. Taking into account the uncertainty in distance as a distribution of distances, 
p{r), it can be shown that the total likelihood in the luminosity domain is 

p{\0g 5*4, n, Sohs\N, fi, (T, Synin, r) = p(l0g Si,n, 5obs|^^, A*S, CrS, Snun), (2.6) 

where jj, and ns are related additively by the term 21ogr, and a and as are equal. The 
final joint posterior in luminosity is then given by 

p{N,fi, cr, S'min, ''I log Si,n, 5obs) ^2 7) 

cx p{logSt,n,Sohs\N,fi,a,S,nin,r) p{N) p{n) p{a) p{Snun) p{r)- 

The prior on A'^ is taken to be uniform from n to cx). We also use uniform priors on the 
model parameters /j, and a. We choose a uniform prior on Sy^in in the range (0, min(S'i)], 
where the upper limit is the flux density of the least bright pulsar in the cluster. The 
prior on r is taken to be a Gaussian. This joint posterior is integrated over various sets 
of model parameters to obtain marginalized posteriors. 



3. Applications 

We applied our Bayesian technique to Terzan 5. Although Terzan 5 has 34 known 
pulsars (Ransom S. M., private communication), we take n = 25, the number of pulsars 
for which we have flux density measurements. The flux densities of the individual pulsars 
were collected in a literature survey (Bagchi et al. 2011 and references therein). The 
flux densities we used were scaled from those reported at 1950 MHz by Ransom et al. 
(2005) and Hessels et al. (2006) to 1400 MHz using a spectral index, a = —1.9, using 
the power law S{i') cx j/". The observed diffuse flux density at 1400 MHz is taken to be 
<5obs = 5.2 mJy (Fruchter & Goss 2000). The prior on N was chosen to be uniform in 
[n, 500], which is sufflciently wide to ensure that the posterior does not rail against the 
prior boundaries. We chose uniform distributions in the same range of 11 and a as used 
by Bagchi et al. (2011) as our priors. We took 5'min to be uniform in (0, min(S'i)]. The 
most recent measurement of the distance to Terzan 5, r — 5.5 ± 0.9 kpc (Ortolani et al. 
2007), was used to model the distance prior as a Gaussian. Figure [T] shows the results of 
the analysis. The median values of the three parameters with 95.45% credible intervals 
are: N = 92+^f , = -0.9+};? and a = 0.9^°;^. 

Note that N is the size of the population of pulsars that are beaming towards the Earth. 
Uncertainties notwithstanding, the beaming fraction of MSPs is generally thought to be 
> 50% (Kramer et al. 1998). This, together with the fact that most pulsars in globular 
clusters are MSPs, imply that the true population size in a cluster is approximately a 
factor of two more than the potentially observable population size. 

3.1. Using prior information 

In the framework developed in the previous section, wc use broad uniform (non-informative) 
priors for the mean and standard deviation of the log-normal. This lack of prior informa- 
tion is apparent in Figure [TJb), where N is not very well constrained. Prior information 
can help better constrain the parameters of interest. Boyles et al. (2011) use models 
of Galactic pulsars from Ridley & Lorimer (2010) to narrow down /i to between —1.19 
and —1.04, and cr to the range 0.91 to 0.98. We have chosen our priors on and cr 
to be uniform within these ranges. Applying the Bayesian analysis over this narrower 
range of /i and a results in much tighter constraints on N as seen in Figure [T]^c), where 
N = 133^53^ This result is consistent with that of Bagchi et al. (2011). 
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Figure 1. Results of the analysis for Terzan 5. For (a) and (b), the analysis was run with wide 
priors on fi and a, with the ranges equal to those used by Bagchi et al. (2011), their Figure 2. (a) 
depicts the joint posterior on /i and a, marginalized over A*', ^min and r; (b) is the marginalized 
posterior for A^. (c) Posterior on after applying narrow priors on /x and a. In (b) and (c), the 
shaded regions lie outside a 95.45% credible interval. 



4. Conclusions 

The technique described here would be useful in future studies of the globular cluster 
luminosity function where ongoing and future pulsar surveys are expected to provide a 
substantial increase in the number of known pulsars in many clusters. We anticipate that 
the increased amount of data would enable us to constrain the distributions of ^ and 
a independently (i.e. without the need to assume prior information from the Galactic 
pulsar population). Further interferometric measurements of the diffuse radio flux in 
many clusters could provide improved constraints on fi and a by measuring the flux 
contribution from the individually unresolvable population of pulsars. 
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